Loci associated with cave-derived traits concentrate in specific regions of the Mexican cavefish genome.

A major goal of modern evolutionary biology is connecting phenotypic evolution with its underlying genetic basis. The Mexican cavefish (Astyanax mexicanus), a characin fish species comprised of a surface ecotype and a cave-derived ecotype, is well suited as a model to study the genetic mechanisms underlying adaptation to extreme environments. Here we map 206 previously published quantitative trait loci (QTL) for cave-derived traits in A. mexicanus to the newest version of the surface fish genome assembly, AstMex3. This analysis revealed that QTL cluster in the genome more than expected by chance, and this clustering is not explained by the distribution of genes in the genome. To investigate whether certain characteristics of the genome facilitate phenotypic evolution, we tested whether genomic characteristics, such as highly mutagenic CpG sites, are reliable predictors of the sites of trait evolution but did not find any significant trends. Finally, we combined the QTL map with previously collected expression and selection data to identify a list of 36 candidate genes that may underlie the repeated evolution of cave phenotypes, including rgrb which is predicted to be involved in phototransduction. We found this gene has disrupted exons in all non-hybrid cave populations but intact reading frames in surface fish. Overall, our results suggest specific "evolutionary hotspots" in the genome may play significant roles in driving adaptation to the cave environment in Astyanax mexicanus and demonstrate how this compiled dataset can facilitate our understanding of the genetic basis of repeated evolution in the Mexican cavefish.

Phylogeographic relationships and morphological evolution between cave and surface Astyanax mexicanus populations (De Filippi 1853) (Actinopterygii, Characidae).

The Astyanax mexicanus complex includes two different morphs, a surface- and a cave-adapted ecotype, found at three mountain ranges in Northeastern Mexico: Sierra de El Abra, Sierra de Guatemala and Sierra de la Colmena (Micos). Since their discovery, multiple studies have attempted to characterize the timing and the number of events that gave rise to the evolution of these cave-adapted ecotypes. Here, using RADseq and genome-wide sequencing, we assessed the phylogenetic relationships, genetic structure and gene flow events between the cave and surface Astyanax mexicanus populations, to estimate the tempo and mode of evolution of the cave-adapted ecotypes. We also evaluated the body shape evolution across different cave lineages using geometric morphometrics to examine the role of phylogenetic signal versus environmental pressures. We found strong evidence of parallel evolution of cave-adapted ecotypes derived from two separate lineages of surface fish and hypothesize that there may be up to four independent invasions of caves from surface fish. Moreover, a strong congruence between the genetic structure and geographic distribution was observed across the cave populations, with the Sierra de Guatemala the region exhibiting most genetic drift among the cave populations analysed. Interestingly, we found no evidence of phylogenetic signal in body shape evolution, but we found support for parallel evolution in body shape across independent cave lineages, with cavefish from the Sierra de El Abra reflecting the most divergent morphology relative to surface and other cavefish populations.

Genetic mapping of craniofacial traits in the Mexican tetra reveals loci associated with bite differences between cave and surface fish.

BACKGROUND: The Mexican tetra, Astyanax mexicanus, includes interfertile surface-dwelling and cave-dwelling morphs, enabling powerful studies aimed at uncovering genes involved in the evolution of cave-associated traits. Compared to surface fish, cavefish harbor several extreme traits within their skull, such as a protruding lower jaw, a wider gape, and an increase in tooth number. These features are highly variable between individual cavefish and even across different cavefish populations. RESULTS: To investigate these traits, we created a novel feeding behavior assay wherein bite impressions could be obtained. We determined that fish with an underbite leave larger bite impressions with an increase in the number of tooth marks. Capitalizing on the ability to produce hybrids from surface and cavefish crosses, we investigated genes underlying these segregating orofacial traits by performing Quantitative Trait Loci (QTL) analysis with F2 hybrids. We discovered significant QTL for bite (underbite vs. overbite) that mapped to a single region of the Astyanax genome. Within this genomic region, multiple genes exhibit coding region mutations, some with known roles in bone development. Further, we determined that there is evidence that this genomic region is under natural selection. CONCLUSIONS: This work highlights cavefish as a valuable genetic model for orofacial patterning and will provide insight into the genetic regulators of jaw and tooth development.

Developmental environment contributes to rapid trait shifts among newly colonized subterranean habitats.

Recent colonization of extreme environments provides unique opportunities to study the early steps of adaptation and the potential for rapid convergent evolution. However, phenotypic shifts during recent colonization may also be due to plasticity in response to changes in the rearing environment. Here, we analyzed a suite of morphological and behavioral traits in paired surface, subterranean, and facultatively subterranean Mexican tetras ( Astyanax mexicanus) from recent introductions in two separate watersheds outside of their native range. We found a variety of phenotypic and behavioral shifts between subterranean and surface populations that are similar to those observed in relatively ancient populations in Mexico. Despite this rapid morphological divergence, we found that most of these trait differences were due to plasticity in response to rearing environments. While most trait assays in common-garden, lab-raised fish indicated that phenotypic shifts in wild fish were the result of plasticity, we also found evidence of genetic control in several traits present in subterranean populations. Interestingly, wall-following behavior, an important subterranean foraging behavior, was greater in lab-born subterranean fish than in lab-born surface fish, suggesting rapid divergence of this trait between subterranean and surface populations. Thus, this study sheds light on the early steps of subterranean evolution, identifies potential rapid behavioral evolution, and suggests that plasticity in traits involving exploratory behavior may facilitate subterranean colonization.

Selection-driven trait loss in independently evolved cavefish populations.

Laboratory studies have demonstrated that a single phenotype can be produced by many different genotypes; however, in natural systems, it is frequently found that phenotypic convergence is due to parallel genetic changes. This suggests a substantial role for constraint and determinism in evolution and indicates that certain mutations are more likely to contribute to phenotypic evolution. Here we use whole genome resequencing in the Mexican tetra, Astyanax mexicanus, to investigate how selection has shaped the repeated evolution of both trait loss and enhancement across independent cavefish lineages. We show that selection on standing genetic variation and de novo mutations both contribute substantially to repeated adaptation. Our findings provide empirical support for the hypothesis that genes with larger mutational targets are more likely to be the substrate of repeated evolution and indicate that features of the cave environment may impact the rate at which mutations occur.

Gene loss and relaxed selection of plaat1 in vertebrates adapted to low-light environments.

Gene loss is an important mechanism for evolution in low-light or cave environments where visual adaptations often involve a reduction or loss of eyesight. The plaat gene family are phospholipases essential for the degradation of organelles in the lens of the eye. They translocate to damaged organelle membranes, inducing them to rupture. This rupture is required for lens transparency and is essential for developing a functioning eye. Plaat3 is thought to be responsible for this role in mammals, while plaat1 is thought to be responsible in other vertebrates. We used a macroevolutionary approach and comparative genomics to examine the origin, loss, synteny, and selection of plaat1 across bony fishes and tetrapods. We show that plaat1 (likely ancestral to all bony fish + tetrapods) has been lost in squamates and is significantly degraded in lineages of low-visual acuity and blind mammals and fish. Our findings suggest that plaat1 is important for visual acuity across bony vertebrates, and that its loss through relaxed selection and pseudogenization may have played a role in the repeated evolution of visual systems in low-light-environments. Our study sheds light on the importance of gene-loss in trait evolution and provides insights into the mechanisms underlying visual acuity in low-light environments.

Widespread convergent evolution of alpha-neurotoxin resistance in African mammals.

Convergent evolution is central to the study of adaptation and has been used to understand both the limits of evolution and the diverse patterns and processes which result in adaptive change. Resistance to snake venom alpha-neurotoxins (αNTXs) is a case of widespread convergence having evolved several times in snakes, lizards and mammals. Despite extreme toxicity of αNTXs, substitutions in its target, the nicotinic acetylcholine receptor (nAChR), prevent αNTX binding and render species resistant. Recently, the published meerkat (Herpestidae) genome revealed that meerkats have the same substitutions in nAChR as the venom-resistant Egyptian mongoose (Herpestidae), suggesting that venom-resistant nAChRs may be ancestral to Herpestids. Like the mongoose, many other species of feliform carnivores prey on venomous snakes, though their venom resistance has never been explored. To evaluate the prevalence and ancestry of αNTX resistance in mammals, we generate a dataset of mammalian nAChR using museum specimens and public datasets. We find five instances of convergent evolution within feliform carnivores, and an additional eight instances across all mammals sampled. Tests of selection show that these substitutions are evolving under positive selection. Repeated convergence suggests that this adaptation played an important role in the evolution of mammalian physiology and potentially venom evolution.

Assessing genomic and ecological differentiation among subspecies of the rough-footed mud turtle, Kinosternon hirtipes.

Combining genetic and ecological measures of differentiation can provide compelling evidence for ecological and genetic divergence among lineages. The rough-footed mud turtle, Kinosternon hirtipes, is distributed from the Trans-Pecos region of Texas to the highlands of Central Mexico and contains 6 described subspecies, 5 of which are extant. We use ddRAD sequencing and species distribution models to assess levels of ecological and genetic differentiation among these subspecies. We also predict changes in climatically suitable habitat under different climate change scenarios and assess levels of genetic diversity and inbreeding within each lineage. Our results show that there is strong genetic and ecological differentiation among multiple lineages within K. hirtipes, and that this differentiation appears to be the result of vicariance associated with the Trans-Mexican Volcanic Belt. We propose changes to subspecies designations to more accurately reflect the evolutionary relationships among populations and assess threats to each subspecies.

CaveCrawler: an interactive analysis suite for cavefish bioinformatics.

The growing use of genomics in diverse organisms provides the basis for identifying genomic and transcriptional differences across species and experimental conditions. Databases containing genomic and functional data have played critical roles in the development of numerous genetic models but most emerging models lack such databases. The Mexican tetra, Astyanax mexicanus exists as 2 morphs: surface-dwelling and cave-dwelling. There exist at least 30 cave populations, providing a system to study convergent evolution. We have generated a web-based analysis suite that integrates datasets from different studies to identify how gene transcription and genetic markers of selection differ between populations and across experimental contexts. Results of diverse studies can be analyzed in conjunction with other genetic data (e.g. Gene Ontology information), to enable biological inference from cross-study patterns and identify future avenues of research. Furthermore, the framework that we have built for A. mexicanus can be adapted for other emerging model systems.

Hybridization underlies localized trait evolution in cavefish.

Introgressive hybridization may play an integral role in local adaptation and speciation (Taylor and Larson, 2019). In the Mexican tetra Astyanax mexicanus, cave populations have repeatedly evolved traits including eye loss, sleep loss, and albinism. Of the 30 caves inhabited by A. mexicanus, Chica cave is unique because it contains multiple pools inhabited by putative hybrids between surface and cave populations (Mitchell et al., 1977), providing an opportunity to investigate the impact of hybridization on complex trait evolution. We show that hybridization between cave and surface populations may contribute to localized variation in traits associated with cave evolution, including pigmentation, eye development, and sleep. We also uncover an example of convergent evolution in a circadian clock gene in multiple cavefish lineages and burrowing mammals, suggesting a shared genetic mechanism underlying circadian disruption in subterranean vertebrates. Our results provide insight into the role of hybridization in facilitating phenotypic evolution.

Whole-genome variation of transposable element insertions in a maize diversity panel.

Intact transposable elements (TEs) account for 65% of the maize genome and can impact gene function and regulation. Although TEs comprise the majority of the maize genome and affect important phenotypes, genome-wide patterns of TE polymorphisms in maize have only been studied in a handful of maize genotypes, due to the challenging nature of assessing highly repetitive sequences. We implemented a method to use short-read sequencing data from 509 diverse inbred lines to classify the presence/absence of 445,418 nonredundant TEs that were previously annotated in four genome assemblies including B73, Mo17, PH207, and W22. Different orders of TEs (i.e., LTRs, Helitrons, and TIRs) had different frequency distributions within the population. LTRs with lower LTR similarity were generally more frequent in the population than LTRs with higher LTR similarity, though high-frequency insertions with very high LTR similarity were observed. LTR similarity and frequency estimates of nested elements and the outer elements in which they insert revealed that most nesting events occurred very near the timing of the outer element insertion. TEs within genes were at higher frequency than those that were outside of genes and this is particularly true for those not inserted into introns. Many TE insertional polymorphisms observed in this population were tagged by SNP markers. However, there were also 19.9% of the TE polymorphisms that were not well tagged by SNPs (R2 < 0.5) that potentially represent information that has not been well captured in previous SNP-based marker-trait association studies. This study provides a population scale genome-wide assessment of TE variation in maize and provides valuable insight on variation in TEs in maize and factors that contribute to this variation.

The utility of genomic prediction models in evolutionary genetics.

Variation in complex traits is the result of contributions from many loci of small effect. Based on this principle, genomic prediction methods are used to make predictions of breeding value for an individual using genome-wide molecular markers. In breeding, genomic prediction models have been used in plant and animal breeding for almost two decades to increase rates of genetic improvement and reduce the length of artificial selection experiments. However, evolutionary genomics studies have been slow to incorporate this technique to select individuals for breeding in a conservation context or to learn more about the genetic architecture of traits, the genetic value of missing individuals or microevolution of breeding values. Here, we outline the utility of genomic prediction and provide an overview of the methodology. We highlight opportunities to apply genomic prediction in evolutionary genetics of wild populations and the best practices when using these methods on field-collected phenotypes.

Repeated evolution of circadian clock dysregulation in cavefish populations.

Circadian rhythms are nearly ubiquitous throughout nature, suggesting they are critical for survival in diverse environments. Organisms inhabiting largely arrhythmic environments, such as caves, offer a unique opportunity to study the evolution of circadian rhythms in response to changing ecological pressures. Populations of the Mexican tetra, Astyanax mexicanus, have repeatedly invaded caves from surface rivers, where individuals must contend with perpetual darkness, reduced food availability, and limited fluctuations in daily environmental cues. To investigate the molecular basis for evolved changes in circadian rhythms, we investigated rhythmic transcription across multiple independently-evolved cavefish populations. Our findings reveal that evolution in a cave environment has led to the repeated disruption of the endogenous biological clock, and its entrainment by light. The circadian transcriptome shows widespread reductions and losses of rhythmic transcription and changes to the timing of the activation/repression of core-transcriptional clock. In addition to dysregulation of the core clock, we find that rhythmic transcription of the melatonin regulator aanat2 and melatonin rhythms are disrupted in cavefish under darkness. Mutants of aanat2 and core clock gene rorca disrupt diurnal regulation of sleep in A. mexicanus, phenocopying circadian modulation of sleep and activity phenotypes of cave populations. Together, these findings reveal multiple independent mechanisms for loss of circadian rhythms in cavefish populations and provide a platform for studying how evolved changes in the biological clock can contribute to variation in sleep and circadian behavior.

Pleiotropic function of the oca2 gene underlies the evolution of sleep loss and albinism in cavefish.

Adaptation to novel environments often involves the evolution of multiple morphological, physiological, and behavioral traits. One striking example of multi-trait evolution is the suite of traits that has evolved repeatedly in cave animals, including regression of eyes, loss of pigmentation, and enhancement of non-visual sensory systems.1,2 The Mexican tetra, Astyanax mexicanus, consists of fish that inhabit at least 30 caves in Mexico and ancestral-like surface fish that inhabit the rivers of Mexico and southern Texas.3 Cave A. mexicanus are interfertile with surface fish and have evolved a number of traits, including reduced pigmentation, eye loss, and alterations to behavior.4-6 To define relationships between different cave-evolved traits, we phenotyped 208 surface-cave F2 hybrid fish for numerous morphological and behavioral traits. We found differences in sleep between pigmented and albino hybrid fish, raising the possibility that these traits share a genetic basis. In cavefish and other species, mutations in oculocutaneous albinism 2 (oca2) cause albinism.7-12 Surface fish with mutations in oca2 displayed both albinism and reduced sleep. Further, this mutation in oca2 fails to complement sleep loss when surface fish harboring this engineered mutation are crossed to independently evolved populations of albino cavefish with naturally occurring mutations in oca2. Analysis of the oca2 locus in wild-caught cave and surface fish suggests that oca2 is under positive selection in 3 cave populations. Taken together, these findings identify oca2 as a novel regulator of sleep and suggest that a pleiotropic function of oca2 underlies the adaptive evolution of albinism and sleep loss.

A chromosome-level genome of Astyanax mexicanus surface fish for comparing population-specific genetic differences contributing to trait evolution.

Identifying the genetic factors that underlie complex traits is central to understanding the mechanistic underpinnings of evolution. Cave-dwelling Astyanax mexicanus populations are well adapted to subterranean life and many populations appear to have evolved troglomorphic traits independently, while the surface-dwelling populations can be used as a proxy for the ancestral form. Here we present a high-resolution, chromosome-level surface fish genome, enabling the first genome-wide comparison between surface fish and cavefish populations. Using this resource, we performed quantitative trait locus (QTL) mapping analyses and found new candidate genes for eye loss such as dusp26. We used CRISPR gene editing in A. mexicanus to confirm the essential role of a gene within an eye size QTL, rx3, in eye formation. We also generated the first genome-wide evaluation of deletion variability across cavefish populations to gain insight into this potential source of cave adaptation. The surface fish genome reference now provides a more complete resource for comparative, functional and genetic studies of drastic trait differences within a species.

It's time to stop sweeping recombination rate under the genome scan rug.

Different parts of the genome can vary widely in their evolutionary histories and sequence divergence from other species. Indeed, some of the most interesting biology (e.g., hybridization, horizontal gene transfer, variable mutation rates across the genome) is revealed by the discordant relationships between taxa across the genome. The goal for much of evolutionary genetics is centred on understanding the evolutionary processes by which such varied signatures arise and are maintained. Many evolutionary genetics studies seek to identify signatures of positive selection between two closely related ecotypes or taxa by delineating regions with particularly high divergence relative to a genome-wide average, often termed "divergence outliers." In a From the Cover article in this issue of Molecular Ecology, Booker et al. take a major step forward in showing that recombination rate differences are sufficient to create false positive divergence outliers, even under neutrality. They demonstrate that the variance of genome scan metrics is especially high in regions with low recombination rates, consistent with previous work. Furthermore, they show that both relative and absolute measures of divergence (FST and DXY , respectively) as well as other commonly used statistics in genome scans (e.g., πW , Tajima's D and H12) all have similar covariance between variance and local recombination rate. Finally, Booker et al. show that low recombination regions will tend to produce more outliers if genome-wide averages are used as cut-offs to define genomic outliers. Booker et al.'s results suggest that recombination rate variation, even under neutral conditions, can shape genome scans for selection, and this important variable can no longer be ignored.

Dark world rises: The emergence of cavefish as a model for the study of evolution, development, behavior, and disease.

A central question in biology is how naturally occurring genetic variation accounts for morphological and behavioral diversity within a species. The Mexican tetra, Astyanax mexicanus, has been studied for nearly a century as a model for investigating trait evolution. In March of 2019, researchers representing laboratories from around the world met at the Sixth Astyanax International Meeting in Santiago de Querétaro, Mexico. The meeting highlighted the expanding applications of cavefish to investigations of diverse aspects of basic biology, including development, evolution, and disease-based applications. A broad range of integrative approaches are being applied in this system, including the application of state-of-the-art functional genetic assays, brain imaging, and genome sequencing. These advances position cavefish as a model organism for addressing fundamental questions about the genetics and evolution underlying the impressive trait diversity among individual populations within this species.

Unique transcriptional signatures of sleep loss across independently evolved cavefish populations.

Animals respond to sleep loss with compensatory rebound sleep, and this is thought to be critical for the maintenance of physiological homeostasis. Sleep duration varies dramatically across animal species, but it is not known whether evolutionary differences in sleep duration are associated with differences in sleep homeostasis. The Mexican cavefish, Astyanax mexicanus, has emerged as a powerful model for studying the evolution of sleep. While eyed surface populations of A. mexicanus sleep approximately 8 hr each day, multiple blind cavefish populations have converged on sleep patterns that total as little as 2 hr each day, providing the opportunity to examine whether the evolution of sleep loss is accompanied by changes in sleep homeostasis. Here, we examine the behavioral and molecular response to sleep deprivation across four independent populations of A. mexicanus. Our behavioral analysis indicates that surface fish and all three cavefish populations display robust recovery sleep during the day following nighttime sleep deprivation, suggesting sleep homeostasis remains intact in cavefish. We profiled transcriptome-wide changes associated with sleep deprivation in surface fish and cavefish. While the total number of differentially expressed genes was not greater for the surface population, the surface population exhibited the highest number of uniquely differentially expressed genes than any other population. Strikingly, a majority of the differentially expressed genes are unique to individual cave populations, suggesting unique expression responses are exhibited across independently evolved cavefish populations. Together, these findings suggest sleep homeostasis is intact in cavefish despite a dramatic reduction in overall sleep duration.

Using multiple reference genomes to identify and resolve annotation inconsistencies.

BACKGROUND: Advances in sequencing technologies have led to the release of reference genomes and annotations for multiple individuals within more well-studied systems. While each of these new genome assemblies shares significant portions of synteny between each other, the annotated structure of gene models within these regions can differ. Of particular concern are split-gene misannotations, in which a single gene is incorrectly annotated as two distinct genes or two genes are incorrectly annotated as a single gene. These misannotations can have major impacts on functional prediction, estimates of expression, and many downstream analyses. RESULTS: We developed a high-throughput method based on pairwise comparisons of annotations that detect potential split-gene misannotations and quantifies support for whether the genes should be merged into a single gene model. We demonstrated the utility of our method using gene annotations of three reference genomes from maize (B73, PH207, and W22), a difficult system from an annotation perspective due to the size and complexity of the genome. On average, we found several hundred of these potential split-gene misannotations in each pairwise comparison, corresponding to 3-5% of gene models across annotations. To determine which state (i.e. one gene or multiple genes) is biologically supported, we utilized RNAseq data from 10 tissues throughout development along with a novel metric and simulation framework. The methods we have developed require minimal human interaction and can be applied to future assemblies to aid in annotation efforts. CONCLUSIONS: Split-gene misannotations occur at appreciable frequency in maize annotations. We have developed a method to easily identify and correct these misannotations. Importantly, this method is generic in that it can utilize any type of short-read expression data. Failure to account for split-gene misannotations has serious consequences for biological inference, particularly for expression-based analyses.

An Adult Brain Atlas Reveals Broad Neuroanatomical Changes in Independently Evolved Populations of Mexican Cavefish.

A shift in environmental conditions impacts the evolution of complex developmental and behavioral traits. The Mexican cavefish, Astyanax mexicanus, is a powerful model for examining the evolution of development, physiology, and behavior because multiple cavefish populations can be compared to an extant, ancestral-like surface population of the same species. Many behaviors have diverged in cave populations of A. mexicanus, and previous studies have shown that cavefish have a loss of sleep, reduced stress, an absence of social behaviors, and hyperphagia. Despite these findings, surprisingly little is known about the changes in neuroanatomy that underlie these behavioral phenotypes. Here, we use serial sectioning to generate brain atlases of surface fish and three independent cavefish populations. Volumetric reconstruction of serial-sectioned brains confirms convergent evolution on reduced optic tectum volume in all cavefish populations tested. In addition, we quantified volumes of specific neuroanatomical loci within several brain regions that have previously been implicated in behavioral regulation, including the hypothalamus, thalamus, and habenula. These analyses reveal an enlargement of the hypothalamus in all cavefish populations relative to surface fish, as well as subnuclei-specific differences within the thalamus and prethalamus. Taken together, these analyses support the notion that changes in environmental conditions are accompanied by neuroanatomical changes in brain structures associated with behavior. This atlas provides a resource for comparative neuroanatomy of additional brain regions and the opportunity to associate brain anatomy with evolved changes in behavior.